A PROBABILISTIC INTERPRETATION OF A SEQUENCE 
RELATED TO NARAYANA POLYNOMIALS 



TEWODROS AMDEBERHAN, VICTOR H. MOLL, AND CHRISTOPHE VIGNAT 

Abstract. A sequence of coefficients appearing in a recurrence for the Narayana 
polynomials is generalized. The coefficients are given a probabilistic interpreta- 
tion in terms of beta distributed random variables. The recurrence established 
by M. Lasalle is then obtained from a classical convolution identity. Some 
arithmetical properties of the generalized coefficients are also established. 



1. Introduction 

The Narayana polynomials 

r 

(1.1) A/;(z) = ^iV(r,fc)z'=-i 

fe=i 

with the Narayana numbers N{r, k) given by 

have a large number of combinatorial properties. In a recent paper, M. Lasalle |19j 
established the recurrence 

(1.3) (Z + l)Mr{z) - Mr+l{z) = Y,{-zY(^^^^\A,,Nr-2n+l{z). 

n>l ^ ^ ' 

The numbers An satisfies the recurrence 

Tl-l 



(1-4) (-1)"-^A„ = C„ + ^(-1)^ iV^-^"-^' 



with A\ = \ and C„ = (^"^ the Catalan number. This recurrence is taken here 

^ n+l V n / 

as being the definition of An ■ The first few values are 

(1.5) Ai = 1, ^2 = 1, As = 5, A4 = 56, A5 ^ 1092, Aq = 32670. 

Lasalle [19] shows that {An : n € N} is an increasing sequence of positive 
integers. In the process of establishing the positivity of this sequence, he contacted 
D. Zeilberger, who suggested the study of the related sequence 

(1-6) an = -p;-, 
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j - 1/ \j + IJ n - j + 1 



with first few values 

(1.7) ai = 2, 02 = 1, as = 2, 04 = 8, 05 = 52, ae = 495, ay = 6470. 
The recurrence ()1.4p yields 

(1.8) (-l)"-ia„ = 2 + ^(-l--^" '^^"+' 

This may be expressed in terms of the numbers 
, , 2 / n \ fn + l 

(1-9) ^n,r--^-[ J 

rt \r — 1/ \r + 1 

that appear as entry ^108838 in OEIS and count Dyck paths by the number of 
long interior inclines. The fact that t7„ ^ is an integer also follows from 

-(::2)(::: 

The relation (11.81) can also be written as 



-l n— 1 

2 + -Y,{-iya^.,c 



(1.11) a„-(-l) 

^ i=i 

The original approach by M. Lasalle [12] is to establish the relation 

(1.12) {Z + \)K{Z) ^ Ur+l{z) = Y,{-zT(^^^\ArXr)Mr^2n+l{z) 

for some coefhcient A„(r). The expression 

(1.13) Nr{z) - ^ Z"(^ + l)'-2™-l ('''^"^^') 

m>0 ^ ^ 

given in [H], is then employed to show that A„(r) is independent of r. This is the 
definition of An given in [Tnj. Lasalle mentions in passing that "J. Novak observed, 
as empirical evidence, that the integers (— 1)"~^A„ are precisely the (classical) 
cumulants of a standard semicircular random variable" . 

The goal of this paper is to revisit Lasalle's results, provide probabilistic inter- 
pretation of the numbers An and to consider Zeilberger's suggestion. 

The probabilistic interpretation of the numbers An starts with the semicircular 
distribution 



fVl -a;2 if - 1 < X < 1 
otherwise. 

Let X be a random variable with distribution /i. Then X* = 2X satisfies 

(1.15) EK] = |0 

I Cm. II r IS even, with r = Zm, 

where C„ = — W (^™) are the Catalan numbers. The moment generating function 

m+l \ m / GO 



(1.14) fi{x) = 



^ 4-n 

(1.16) ^{t)^Y.^[X-]- 

n=0 
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is expressed in terms of the modified Bessel function of tlie first kind (cc) and tlie 
cumulant generating function 

(1.17) V(i)=log¥'(i) = E^i('^)^ 

n=l 

has coefficients Ki(n), known as the cumulants of X. The identity 

(1.18) A„ = (-1)"+Vi(2n)22", 

is estabUshed here. Lasalle's recurrence (|1.4p now follows from the convolution 
identity 



ri-l 



71-1 



(1.19) nin) = E [X"] ~ [ . [^""^J 

that holds for any pair of moments and cumulants sequences [241. The coefficient 
a„ suggested by D. Zeilberger now takes the form 

(120) a 2(-l)"+^.i(2n) 

In this paper, these notions are extended to the case of random variables dis- 
tributed according to the symmetric beta distribution 

(1.21) f,ix)= ^ ^ a-x'r-'/^ foT\x\<i,f,>~i 

-^IM + 2' 2> 

and otherwise. The semi-circular distribution is the particular case fj, — 1. Here 
B{a, b) is the classical beta function defined by the integral 

(1.22) B{a,b)^ f f-^{l~tf-^ dt, fora,6>0. 

Jo 

These ideas lead to introduce a generalization of the Narayana polynomials and 

these are expressed in terms of the classical Gegenbauer polynomials C„ ^ . The 
coefficients a„ are also generalized to a family of numbers {a„(/^)} with parameter 
^. The special cases — Q and — ±\ are discussed in detail. 

Section[2]produces a recurrence for {an} from which the fact that an is incresaing 
and positive are established. The recurrence comes from a relation between {an} 
and the Bessel function Ia{x). Section [3] gives an expression for {on} in terms of a 
determinant of an upper Hessenberg matrix. The standard procedure to evaluate 
these determinants gives the original recurrence defining {an}. Section[4]introduces 
the probabilistic interpretation of the numbers {a„}. The cumulants of the asso- 
ciated random variable are expressed in terms of the Bessel zeta function. Section 
[5] presents the Narayana polynomials as expected values of a simple function of 
a semicircular random variable. These polynomials are generalized in Section [S] 
and they are expressed in terms of Gegenbauer polynomials. The corresponding 
extension of {an} are presented in Section[7l The paper concludes with some arith- 
metical properties of {an} and its generalization corresponding to the parameter 
^ ~ 0. These are described in Section [H 
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2. The sequence {a„} is positive and increasing 

In this section a direct proof of the positivity of the numbers an defined in (|1.8|) 
is provided. Naturally this implies An > 0. The analysis employs the modified 
Bessel function of the first kind 

Formulas for this function appear in |16| . 
Lemma 2.1. The numbers a„ satisfy 



(2-2) E ■ 



Proof. The statement is equivalent to 

(2.3) Vi/,(2Vi) X E + J-r)l = 2^^(2v/i). 

This is established by comparing coefficients of on both sides and using (|1.8p . □ 
Now change x to in Lemma |2. II to write 

{-ly-^oj .x2J-2 _ 2/2(2a;) 
^tt (.? + !)! ~ ^ 

The classical relations 

(2.5) ^ (z-"/™(z)) = z-^Wilz), and ^ (z™+i/„,+i(z)) - z"+i/™(z) 
give 

(2.6) I[{z)^h{z) + U^{z). 



(2-4) E ■ 



Therefore ()2.4p may be written as 

(2.7) y^MTS^^iAi^gfMM 

^ ^ (j + l)! (j-1)! xdx 2x 

The relations (|2.5p also produce 



V z'^I„,{z) J Il{z) 
In particular, 



Replacing this relation in (|2.7I) gives the recurrence stated next. 
Proposition 2.2. The numbers a„ satisfy the recurrence 

(2.10) 2na„ = E (^^ " (^^ -^^"fc«»-fc' for n > 2, 

with initial condition ai = 1. 
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Corollary 2.3. The numbers a„ are nonnegative. 
Proposition 2.4. The numbers a„ satisfy 

(2.11) 4a„ = J2 (l I I) I ^) - E 

fe=l ^ ^ ^ ^ fe=2 

Proof. This follows from (|2.10p and the identity 
n \ f n \ n /n — l\/n— 1 

fc- 1/ U+ 1/ 2 U - 1 



n — 1 \ / n — 1 
fc-2 



77, — 1 \ /n — 1 



k-2 



□ 



Corollary 2.5. The numbers a„ are nonnegative integers. Moreover a„ is even if 
n is odd. 



Proof. Corollarv 12.31 shows a„ > 0. It remains to show a„ G Z and to verify the 
parity statement. This is achieved by simultaneous induction on n. 

Assume first n ~ 2m + 1 is odd. Then (jl.9|) shows that ^cr„.T. G Z and (|l.lip . 
written as 

n-l 

(2.12) an = (-1)"-^ ' 

proves that a„ £ Z. Now write (12.101) as 
(2.13) 



/2m + l\ /2m + l 
2(2m + l)a2m+i 2 2_^ ( , , ) I , , )afea2m+i-fc 



fe=i 



fc - 1 



fc + 1 



and observe that either k or 2r7i + 1 — A: is odd. The induction hypothesis shows 
that either Ofc or a2m+i-fe is even. This shows a2m+i is even. 

Now consider the case n = 2m even. If r is odd, then ar is even; if r is even then 
r — 1 is odd and ^an,r € Z in view of the identity 

2 /?! - 1\ /n + 1^ 
2)[r+ , 

□ 



(2.14) 



r — 1 \r 

The result follows again from (|2.12p . 
Corollary 2.6. The numbers An are nonnegative integers. 



The recurrence in Proposition 12.21 is now employed to prove that {a„} is an 
increasing sequence. The first few values are 2, 1, 2, 8, 52. 

Theorem 2.7. For n > 3, the inequality a„ > a„_i holds. 

Proof. Take the terms k — \ and fc = n — 1 in the sum appearing in the recurrence 
in Proposition (|2.2p and use a„ > to obtain 



1 

2n 



n \ / ri \ in 



(2.15) , , , , , 

Since ai = 2 the previous inequality yields 

(2.16) a„ > (n - l)a„_i. 
Hence, for n > 3, this gives a„ — a„_i > (n — 2)a„_i > 



□ 
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3. An expression in forms of determinants 
The recursion relation (|1.8|) expressed in the form 



(3.1) 



E 



2m 



is now employed to produce a system of equations for the numbers a„ by varying 
m through 1, 2, 3, • • • ,n. The coefficient matrix has determinant (— l)^^)^! and 
Cramer's rule gives 



(3.2) 



■ dct 



(A) (III) 



\2^lJ \2+l) 



2+17 







\3-lJ \3+lJ 



2\ 

4 

6 



V(i-i)(i:d G-i)(2:d L-,)c;ti) ••• ejrr) w 

The power of —1 is eliminated by permuting the columns to produce the matrix 

\ 



(3.3) B„ 



/2 Gioa:;) o 
6 



V2n (3-1) (aS--- L-2)(T)/ 

The representation of a„ in terms of determinants is given in the next result 

Proposition 3.1. The number a„ is given by 
(3.4) 

where Bn is the matrix in p.3p . 



det B„ 



Recall that an upper Hessenberg matrix is one of the form 
/3i,2 

/32,1 /32,2 /32,3 
P3.1 P3,2 P3.3 P3A 



(3.5) 



Hn — 



\PnS Pn,2 Pn,3 PnA 

The matrix B is of this form with 



\ 






(3.6) 



2i 



(A)Ct') 



if 1 < i < n and j = 1 
iij— l<i<n and j > 1 . 



It turns out that the recurrence (jl.Sp used to define the numbers a„ can be 
recovered if one employs p.4p . 
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Proposition 3.2. Define an by 

„^ detBn 
(3.7) a„ = j — 

where B is the matrix (j3.3l) . Then {a„} satisfies the recursion 

n — 1 \ / n + 1 \ a j 



(3.8) (-l)"-V = 2 + ^(-l)^- , 

and the initial condition ai — 1. Therefore a„ — an 



Proof. For convenience define det = 1. Tlic determinant of a Hessenberg matrix 
satisfies the recurrence 

n n—1 

(3.9) detHn = ^(-l)"-'^/?„,.det77,_i [] 

r— 1 z— r 

A direct application of (|3.9p yields 



1 ( (-l)"-(2n)(n - 1)! + ±(-1)-^ ( ^_ + det 5.., J] 



r=2 ^ ^ ^ 



2(_l)"-i+^(_l)"-' " ^^"+1^ 



r-l 



r — 2/ \ r J n — r + 2 

■+(-i)"-'i:(-i)'(":;)(":;);r^' 



2(-l) 

This is eiHi). □ 



Corollary 3.3. The modified Bessel function of the first kind admits a determinant 
expression 

(-l)J-idctS, /a;\2j 
(3.10) /r(.) = X exp ( g (/+i),^.,2^ (2) 

Proof. This follows by integrating the identity 

□ 



4. The probabilistic background: conjugate random variables 

This section provides the probabilistic tools required for an interpretation of the 
sequence An defined in (|1.4p . The specific connections are given in Section O 
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Consider a random variable X with the symmetric beta distribution given in 
()1.21|) . The moments of the symmetric beta distribution, given by 

1 /■! 

(4.1) E[X"] = — / x'^il-xY'^/^dx, 



B{tJi+\,\) 
vanish for n odd and for n — 2m they are 

r(/i + 1) (2m)! 



(4.2) E [X^™] 



^{^l + l + m) 22™ m!' 
Therefore the moment generating function is 

(4.3) ^Jt) = E [e*^] = y E [X"] — = r(^ + 1) y -. 

n— m— 

The next proposition summarizes properties of ^fi{t). The first one is to recognize 
the series in (|4.3p from (|2.ip . The zeros of the Bessel function of the first 

kind 

(4-4) ,w i ^n (f) 

■■^ m! 1 (m + a + 1) V2/ 

appear in the factorization of ip^^ in view of the relation I^{z) = e^'^^'^/^ J^(iz). 

Proposition 4.1. The moment generating function iffj,{t) of a random variable 
X ^ .fy,is given by 

(4.5) ^,,(t)=r(/i+l) 0^/^(0. 

Note 4.2. The Catalan numbers C„ appear as the even-order moments of when 
ji — 1. More precisely, if X is distributed as /i (written as X ^ /i), then 

(4.6) E [(2X)2"] = C„ and E [(2X)2"+i] = 0. 

Note 4.3. The moment generating function of admits the Weierstrass product 
representation 

fc=i V -^M.fc/ 

where {j^j./t} are the zeros of the Bessel function of the first kind J^. 
Definition 4.4. The cumulant generating function is 

Ip^it) = log if ^{t) 

/ °° j-n\ 

= \og[Y,E[X-] 



\n=0 
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The product representation of [t) yields 
log VpW 



k=i \ -^M.fe / 



oo^ oo^ (_l)n-l / ^ X 2n 



n. 

The series converges for \t\ < j^^i. The first Bessel zero satisfies j^^i > for all 
/i > 0. It follows that the series has a non-zero radius of convergence. 

Note 4.5. The coefficient K^((ri) is the ri-th cumulant of X. An expression that 
links the moments to the cumulants of X is provided by V. P. Leonov and A. N. 
Shiryaev [20]: 

(4.8) n,{n) = ^(-l)^-i(fc - 1)! H E(2X)I^^I 

V i=l 

where the sum is over all partitions V = {Vi, • • • , Vfc} of the set {1, 2, . . . , n}. 

In the case /i = the moments are Catalan numbers or 0, in the case ^ — \ the 
moments are central binomial coefficients. Therefore, in both cases, the cumulants 
K,fi{n) are integers. An expression for the general value of (x involves 

CO 

(4-9) C.(^) = E-^ 

k=i •'tJ'-M 

the Bessel zeta function^ sometimes referred as the Rayleigh function. 

The next result gives an expression for the cumulants of a random variable X 
with a distribution The special case ^ = 1, described in the next section, 
provides the desired probabilistic interpretation of the original sequence An ■ 

Theorem 4.6. Let X - Then 

if ri is odd, 

2(-l)"/2+i(„_ l)!C^(n) if n is even. 

Proof. Rearranging the expansion in Definition 14.41 gives 

2ra 



(4.10) K(n) = 



(-1)"+! / t \ 



n=l fc=l 

Now compare powers of t in this expansion with the definition 
(4.11) log^^(i) = E^^'(")- 

1 

n— 1 

to obtain the result. □ 
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The next ingredient in the search for an interpretation of the sequence An is the 
notion of conjugate random variables. The properties described below appear in 
[25) . A complex- valued random variable Z is called a regular random variable {rrv 
for short) if E|Z|" < oo for all n e N and 

(4.12) E[h{Z)]=h{E[Z]) 

for all polynomials h. The class of rrv is closed under compositions with polynomials 
(if Z is rrv and P is a polynomial, then P{Z) is rrv) and it is also closed under 
addition of independent rrv. The basic definition is stated next. 

Definition 4.7. Let X, Y be real random variables, not necessarily independent. 
The pair {X, Y) is called conjugate random variables if Z — X + iY is an rrv. The 
random variable X is called self- conjugate if Y has the same distribution as X . 

The property of rrv may be expressed in terms of the function 

$(a, 13) E [exp{iaX + if^Y)] 

The next theorem gives a condition for Z = X + iY to he an rrv. The random 
variables X and Y are not necessarily independent. 

Theorem 4.8. Let Z = X+iY be a complex valued random variable with E [Z] = 
and E [Z"] < oo. Then Z is an rrv if and only if $(«, ia) — 1 for all a E C 

This is now reformulated for real and independent random variables. 

Theorem 4.9. Let X, Y be independent real valued random variables with finite 
moments. Define 

$x(a) = E [e'"^] = V ^^E [X"] and $y(/3) = E [e'"^] - V ^^E [F"] . 



Then Z = X + iY is an rrv with mean zero if and only if ^ x {ct)^Y (ict) = 1- 

Example 4.10. Let X and Y be independent Gaussian variables with zero mean 
and the same variance. Then X and Y are conjugate since 

fxit) = exp (^'Y^^^ fivit) = cxp (^^~Y^' 

Note 4.11. Suppose Z ^ X + iY is a rrv with E [Z] = and z eC. The condition 
(|4.12l) becomes 

(4.13) E[h{z + X + iY)] = h{z). 

Given a sequence of polynomials {Qn{z)} such that deg{Qn) ~ n and with 
leading coefficient 1, an elementary argument shows that there is a unique sequence 
of coefficients Uj^n such that the relation 

n 

(4.14) Qn+i{z) - zQn{z) = ^ aj- „Qj(z) 

holds. This section discusses this recurrence for the sequence of polynomials 

(4.15) P„(z) :=E(z + X)" 

associated to a random variable X . The polynomial Pn is of degree n and has 
leading coefficient 1. It is shown that if the cumulants of odd order vanish, then 
the even order cumulants provide the coefficients aj^n for the recurrence (|4.14p . 
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Theorem 4.12. Let X be a random variable with cumulants K{m). Assume the 
odd-order cumulants vanish and that X has a conjugate random variable Y . Define 
the polynomials 

(4.16) P„(z) =E[(z + Xr]. 
Then P„ satisfies the recurrence 

(4.17) Pn+iiz) - zPniz) = > ( , , U(2m)P„_ 

m>l ^ ' 

Proof. Let Xi, X^ independent copies of X. Then 
E [Xi {{Xi +1Y1 + Z + X2T + X2T)] = 



3=0 

This last expression becomes 



^ [xi(Xi + tY,y{z + X2r-^] = E (^)iE [^1(^1 + ^YiV] E [{z + X2r-q 



3 

On the other hand 



E [Xi {{Xi +Z + X2+ ^rl)" -iz + X2)")] = 

J2 f'')lE [Xi(Xi + z)"-^ E [(^2 + ^Yl)1 - E [X,{z + X2)"] . 



The cancellation property (j4.28l) shows that the only surviving term in the sum is 
r = 0, therefore 

E [Xi {{Xi + Z + X2+ ^Yl)" - (z + X2)")] = 

E [Xi{Xi + z)"] - E [Xi] E [{z + X2)"] 
and E [Xi] ~ since k(1) — 0. This shows the identity 

(4.18) (^^'^E[Xi(Xi+zyi)-']E[(z + X2)"-^] =E[Xi(Xi + z)"]. 

The cumulants of X satisfy 

(4.19) K{m) =EX{X + iY)"'-\ for m > 1, 
(see Theorem 3.3 in [13]), therefore in the current situation 

(4.20) E[X,iX,+.Y,y]J' 

Ik(2to) it j = 2m + 1 is odd. 

On the other hand 

E[Xi(Xi+z)"] = E[(Xi+z)"+1-z(Xi+z)"] 

= Pn+l{z) - zPn{z). 

Replacing in ()4.18p yields the result. □ 
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Recall that a random variable has a Laplace distribution if its distribution func- 
tion is 



(4.21) 



Assume has a distribution defined in (|1.2ip and moment generating func- 
tion given by (|4.7p . The next lemma constructs a random variable conjugate to 
X.,.. 



Lemma 4.13. Let y^,„ be a random variable defined by 
(4.22) 



where {Lk : A: S N} is a sequence of independent, identically distributed Laplace 
random variables. Then lim Y,, ^ 



- Kj exists and is a random variable with 

n— f C30 

continuous probability density. Moreover, the moment generating function of lYf^ 
is 

(4.23) E [e^*^-] = n f 1 + -^-) ■ 

fe=l V ^t^.kj 

the reciprocal of the moment generating function of given in (|4.7p . 
Proof. The characteristic function of a Laplace random variable iLk/ji_i,k is 

(4.24) ^ 
The values 

(4.25) E 



1 









= 0, and E 




Jfi,k 



JtJ,,k 



guarantee the convergence of the series 



(4.26) 



OO 

E 

fc=i 



E 



and 



k=l 



3. 

Jfi,k 



(The last series evaluates to 1/(2^ -|- 2)). This ensures the existence of the limit 
defining (see [17] for details) . The continuity of the limiting probability density 
Yfj^ is ensured by the fact that at least one term (in fact all) in the defining sum has 
a continuous probability density that is of bounded variation. □ 

Note 4.14. In the case ~ is independent of F^, then the conjugacy property 
states that if h is an analytic function in a neighborhood O of the origin, then 



(4.27) 

In particular 
(4.28) 



E[h{z + Xf, + iYf,)]^ h{z), ioi-zeO. 



E [{X^ + tY^r] = 



if n = 0, 
otherwise. 



Note 4.15. In the special case — n/2 — 1 for n e N, n > 3, the function (j4.23p 
has been characterized in |llj as the moment generating function of the total time 
Tn spent in the sphere 5""^ by an n-dimensional Brownian motion starting at the 
origin. 
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5. The Narayana polynomials and the sequence An 

The resuh of Theorem 14.121 is now apphed to a random variable X ^ fi. In this 
case the polynomials P„ correspondi, up to a change of variable, to the Narayama 
polynomials Mn- The recurrence established by M. Lasalle comes from the results 
in Section m In particular, this provides an interpretation of the sequence {j4„} in 
terms of cumulants and the Bessel zeta function. 

Recall the distribution function /i 



(5.1) 



2V1 -a;2/7r, 
0, 



for |x| < 1 
otherwise. 



Lemma 5.1. Let X ^ fi. The Narayana polynomials appear as the moments 
(5.2) Afr{z) ^eUi + Z + 2y/^XY' 

for r > 1. 

Proof. The binomial theorem gives 

E [(1 + z + 2^xy-'] = ^ ('^ T + ly-^-^z^/^E [{2xy] 

The result now follows from (|4.6p and (|1.13l) . 

In order to apply Theorem 14. 121 consider the identities 



□ 



(5.3) 



with 
(5.4) 



Ur{z) 



E 



{l + z + 2^X) " 

= {2^y-^e\{x + z.y-^ 

= (2yi)'-ip._i(z.), 



2^ 



The recurrence (|4.17p applied to the polynomial Pn{z^,) yields 



(5.5) 



Un+2{z) {l+z)Mn+l{z) 



E 

m> I 



n 

2m -I 



k(2to) 



Mn- 



-2m+2 



(z) 



(2Vi) 



n— 27n+l 



(2Vi)"+i 2yi {2y^y 
that reduces to 

(5.6) {l + z)K{z)-K+i{z) = -y] (j''^)ti{2m)2^'^z^K+i-2,n{z), 
by using r = n + 1 . This recurrence has the form of (|1.12l) . 

Theorem 5.2. Let X ~ /i. Then the coefficients An in Definition 11.31 are given 

by 

(5.7) An = (-l)"+iK(2n)22". 
The expression in (|4.10p gives the next result. 



14 T. AMDEBERHAN, V. MOLL, AND C. VIGNAT 

Corollary 5.3. Let 

CO 

(5.8) C.(«) = E-^ 

fc=l "'M,fe 

be the Bessel zeta function. Then the coefficients An are given by 

^„ = 22"+i(2n-l)!Ci(2n). 
The scaled coefficients a„ are now expressed in terms of the Bessel zeta function. 
Corollary 5.4. The coefficients a„ are given by 

(5.9) an - 22"+i(7i+ f)!(n- f)!Ci(27i). 
Note 5.5. This expression for the coefficients and the recurrence 

n-l 

(5.10) (n + fi)U2n) = ^ U2r)U2n - 2r). 

r=l 

given in |15| . provides a new proof of the recurrence in Proposition (j2.2p . 

6. The generalized Narayana polynomials 

The Narayama polynomials Afr{z), defined in (|l.ip . have been expressed as the 
moments 

(6.1) K{z) + Z + 2y/^XY' 

for r > 1. Here X is a random variable with distribution function /i. This suggests 
the extension 

(6.2) AA/;(2) =E [(l + z + 2V^X)"" 
with X - /p. Therefore, 7V„ = 7V„i. 

Note 6.1. The same argument given in (j5.6p gives the recurrence 

(6.3) (1 + z)AA/'(z) - ACV,(z) = -J2 (2^„l\)'^(2H2''"^"ACVi-2™(^), 

m>l ^ ' 

where Ac(2n) are the cumulants of X ^ Theorem 15 . 2 1 gives an expression for the 
generalization of the Lasalle numbers: 

(6.4) (-l)"+iK(2n)22" 

and the corresponding expression in terms of the Bessel zeta function: 

(6.5) Aj::=22"+i(2n-l)!C^(2n). 

The generalized Narayana polynomials are now expressed in terms of the Gegen- 
bauer polynomials C!^[x) defined by the generating function 

oo 

(6.6) "^C^{x)f^ ^{\-2xt^eY^. 

n=0 
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These polynomial admit several hypergeometric representations: 



3.7) C^(.t) = 



2"(M)r 

Til 



2F1 —n, n + 2fi; fi + ^ 



1-x 



{x - 1)" 2F1 i-n,-n-^i+ i; -2?! -2^ + 1; ^— 



The connection between Narayana and Gegenbauer polynomials comes from the 
expression for C!^{z) given in the next proposition. 



Proposition 6.2. The Gegenbauer polynomials are given by 



(6.8) 



I (cos6' + isin6'cos(/))"sin^^"^ 
Jo 



Proof. The Laplace integral representation 

appears as Theorem 6.7.4 in P]. The change of variables z = cos0 and X = cost 
gives 

r(n + 2^) 



22'^n!r2(^) 

(Mile 



y" (z + - 1^)" (1 - x^y-^ dx 



as claimed. Since this is a polynomial identity in z, it can be extended to all 

z e C. □ 

Theorem 6.3. The Gegenbauer polynomial C}^ and the generalized polynomial 
Af^ satisfy the relation 



(6.10) 



+ l)n 



Proof. Introduce the variable 

(6.11) 

so that 

(6.12) ; 
Then 



Z - 1 



Z 



and 



(1 - zrc: 

l + z 
' l-z 

z 



+ i f l + z 
1 - z 



l + z 



Z+l ^/Z^^ 2Vi' 



f 1 + z 
1 - z 



(2^ + l)„ / 2^1 



1-z 



E 



( 



i^^E[(l + . + 2ViX,)"] 



n\ (1 - z) 

(2/^ + 1) 
n\ (1 - z)'' 



using - 1 = 4z/(l - z)^. 



□ 
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The expression (|6.7p now provides hypergeometric expressions for the original 
Narayana polynomials 

(6.13) AA„+,(z) = , 'i''"^" (\^) ■ 

{n + 2j{n+l) \1 - z J 

Corollary 6.4. The Narayana polynomials are given by 

(6.14) A/^+iCz) = (i_z)"2Fi (^-n,n + 3;2;-^^ 

(2n + 2)! „ / z-1 



(n + 2)!(n+ 1)! 
= 2^'l(-?^, 1;2;2:). 

This yields the representation as finite sums 



fc=0 

„ , 1 „ , 1 



+ U + 1, 



n + l\ /n + 1 
fc 

Note that the first two expressions coincide up to the change of summation 
variable k —> n — k while the third identity is nothing but (|l.ip . 

Note 6.5. The representation 

(6.16) Ciiiz) = ^-^{2xr,F, (^-^, i^; 1 - „ - ^; ^ 

that appears in as 6.4.12 in Sj, gives the expression 
(6.17) 

, (2n + 2)! /l + z\" / n 1-n 1 /1-z^^^ 

•^"+^(")= (n + l)!(n + 2)! ) l^" 2 ' ~" ^ 2' (,TT^ 

equal to the finite sum representation 



n-2fe 



Note 6.6. The polynomials S'n(z) = zA/';[(z) satisfy the symmetry identity 

(6.19) Sn{z)^Z^+^Sniz~^). 

These polynomials were expressed in [21] as 

(6.20) Sn{z) = (z - 1)"+^ / Pn{2x - 1) dx 



1/2 

where Pn{x) — Cn [x) are the Legendre polynomials. An equivalent formulation 
is provided next. 
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Theorem 6.7. The polynomials Sn{z) are given by 
Proof. The integration rule 

(6.21) J C^{x) dx = -^l—^C':-l{:c) 
implies 

(6.22) f^'^^'~'' C]l\2x -l)dx^ -\C-1{' ( i±j) , 
since the generating function 

oo 

(6.23) ^rC~i/2(2) = (l-2zi + i2)i/2 

gives C^+(^(-l) ^ for n>l. Then (jg:^ yields 

(6.24) 5„(.) = -i(z-l)"+iC-^f 

A classical formula for the Gegenbauer polynomials states 

Ln/2J , / s 

(6.25) C„nz) = 5: falTfcV ^-)"'" 
and the identity 

produce 



n+l-2fc 



1\ 1 (2fc-2)! 



2)^ 22*^-1 (fc-1)! 

L"/2J 



fc=0 ^ / \ / 

The result now follows from (|6.24p . □ 

7. The generalization of the numbers a„ 
The terms forming the original suggestion of Zeilberger 

(7.1) a„ 



Cn 

have been given a probabilistic interpretation: let X be a random variable with 
a symmetric beta distribution function with parameter fJ. = 1 given explicitly in 
(15. ip . The numerator An is 

(7.2) An = {-l)"+^K{2n)2^" 

where K(2n) is the even-order cumulant of the scaled random variable = 2X. 
The denominator C„ is interpreted as the even-order moment of X^, : 

(7.3) Cn=E [^2"] . 

These notions are used now to define an extension of the coefficients a„. 



18 



T. AMDEBERHAN, V. MOLL, AND C. VIGNAT 



Definition 7.1. Let X be a random variable with vanishing odd cumulants. The 
numbers a„(/i) are defined by 

/ N 2(-l)"+iK(2ri) 

(7.4) = 

In the special case X^, — 2X with X ^ f^, these numbers are computed using the 
cumulants 

(7.5) «^(2n) = (-l)"+i22"+i(2n - iy.U2n) 
and the even order moments 

to produce 

(7.7) a„(Ai) = 22"+i {n - 1)! + 1)„ Cp(2n). 

The value 

= ioITT) 

yields the initial condition ai(/i) = 2. 

The recurrence (|5.10p now provides the next result. Recall that when x is not 
necessarily a positive integer, the binomial coefficient is given by 

(7 9) h = ^(" + ^) 

^ ' ' \kj T{x-k + l)kl' 

Proposition 7.2. The coefficients a„(/i) satisfy the recurrence 

(7.10) Q»(a^) = ^^n+M-i) 2^ (^,, _ _ 1 j jafe(M)a„-fe(M), 

with initial condition ai(/i) = 2. 

Proof. Start with the convolution identity for Bessel zeta functions (|5.10p and re- 
place each zeta function by its expression in terms of a„(/i) from (j7.7p . which gives 



l" + M^22"+i(n-l)!(M + l)„ 



ak{^J.) an-fe(M) 



E 22fe+i(fc _ + 1)^ 22»-2'c+i(„ _ A: - 1)!(m + l)„-fc 



and after simplification 



o-nW = TTi — \ — V Ti TTT7 1 — TTT7 — TTTl — rT\ ak(p-)an^k[fJ-). 

2(n + m) ^ (fc - !)'(" -k-iy.{fi + l)k{fJ- + l)n-k 

The resut now follows by elementary algebra. □ 

Note 7.3. In the case /i = 1, the recurrence (|7.10p becomes (|2.10l) and the coeffi- 
cients a„(l) are the original numbers a„. 
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Note 7.4. The recurrence (|7.10|) can be written as 

""('^^ = 2 ^ r(^ + fc + i)r(n + M - fc + mn - k)m «'^^(^)«"-^-('^)- 

Theorem 7.5. The coefBcients a„{fi) are positive and increasing for n > [^^J. 



Proof. The positivity is clear from ()7.7p . Now take the terms corresponding to 
k = 1 and A; = n — 1 in (|7.10p to obtain 

n — 1 2{n — 1) 

(7.11) an(/i) > — — -ai(^)a„_i(/i) = — — a„_i(^). 

This yields 

(7.12) a„(^) - a„_i(^) > ^ ^ a„-i(/x) 

and the result follows. □ 

Some other special cases are considered next. 

The case /i = 0. In this situation the distribution is the arcsine distribution given 
by 

(7.13) /oW = l^^' '"'-'l^' 

10, otherwise. 

By the recurrence on the Co function, the coefficients 

(7.14) a„(0) = 22"(n-l)!n!Co(2n) 
satisfy the recurrence 

(7.15) a„(0) = i J2 fc ^) (l I J)afc(0)«n-;o(0) 

with ai(0) = 2. Now define as Lasalle 6„ = ^an{0) and then (|7.15p becomes 

(7.0) K ^ 

6i = 1. 
In particular 6„ is a positive integer. 

The following comments are obtained by an analysis similar to that for a„. 
Note 7.6. The recurrence 



gives the generating function 



^(-ir^ x^-^ _ /i(2x) _ 1 d 
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Note 7.7. The sequence bn adniits a, determiiiaiit representation bn — det(Af7i), 
where 



(7.17) 



M„. 



1 (?)(?:!) 







\ 



Vi QiTA) QiV,) ••• L-i)C:2)/ 

Note 7.8. The identity 12(2;) = loix) — is expressed as 



(7.18) 

provides the relation 
(7.19) 



xh{2x) 



1 + T^a; 



2/2(23;) 
2"" xli{2x) 



= 1 



J- 



The case = 5- In this situation the distribution is the uniform distribution on 
[— 1, 1] with even moments 

22n 

(7.20) EX^" = 

^ ' * 2n+l 

and vanishing odd moments. The sequence of cumulants is 

(7.21) A«i/2(2n) = 2(-l)»+i(2n - 1)! Ci/2(2n) 
where the Bessel zeta function is 



(7.22) 



°° 1 1 2^""^ 

Ci/2(2n) = = -^C(2n) = ^|S2„| 



(2n)! 



where B„ are the Bernoulli numbers. This follows from the identity 

(7.23) 
This yields 
(7.24) 



Ji/2{x) = \l — sinx. 



K;i/2(2n) = 2^"^ and /ti/2(2n + 1) = 0, 



with Ki/2{0) = 0. These are the coefficients of u"'/n\ in the cumulant moment 
generating function 



(7.25) 



l0gVl/2(w) = log 



sinhu 



1 

I80' 



u 6 

Finally, the corresponding sequence 

2(-l)"+iK(2n) 



2835 



(7.26) 
is given by 
(7.27) 



an (^) = 



E[X2"] 



a„(|)=2^"^|B2„| 
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(7.31) /_i/2(x) = -J(x-l) + -<5(x + l) 



The first few terms are 

(7.28) ai (i) = 2, a2 (i) = | (i) = ^, a4 (i) = as (i) = 

as expected, this is an increasing sequence for n > 3. The convolution identity 
(|5.10p for Bessel zeta functions gives the well-known quadratic relation for the 
BernoulH numbers 

"^1 /2n\ 

(7.29) H\2k] ^2fe^2«-2fe = -(2n + 1)S2„, for n > 1. 

Moreover, the moment-cumulants relation (|1.19p gives, replacing n by 2n and after 
simplification, the other well-known identity 

(7.30) (^^'^l ^) 2''^2j - 2n, for n > 1. 

Note 7.9. The generating function of the sequence a„ (i) is given by 
^3/2(3^) _ xtanhx - 1 _ (-l)J~^2aj- {\) 

The limiting case /it = — 5 has the probability distribution 

(the discrete Rademachcr distribution). For a Rademacher random variable X, the 
odd moments of X* — 2X vanish while the even order moments are 

(7.32) E [X^"] = 2^". 
Therefore 

(7.33) «-i/2(2n) = (-l)"+42"+i(2n - 1)! C-i/2(2n). 
The identity 

(7.34) J_i/2(x) = W — coss 

V T^X 

shows that jk, -1/2 — (2fc — l)''''/2 and therefore 

°° 92" 2^" — 1 

(7-35) C-i/2(2n) = g ^2„(2fc,i)2.„ = ^C(2.). 

The expression for k_2/2(2't-) may be simplified by the relation 

(7.36) ^„ ^ ?^(2"+' - l)S„+i 

n -|- 1 

between the Euler numbers En and the Bernoulli numbers. It follows that 

(7.37) K_i/2(2ri) = -24"-iS2„-i. 
The corresponding sequence a„ (—5) is now given by 

(7.38) a„ (-i) = (-l)"22"£;2„_i 
and its first few values are 

ai (-i) = 2, 02 (-^) = 4, 03 (-5) = 32, a4 (-i) - 544, 05 (-5) = 15872, 
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xl-i/2{x) X j^^ (2j-l)! 



Note 7.11. The convolution identity (jS.lOp yields the well-known quadratic recur- 
rence relation 



and the momcnt-cumulant relation (|1.19p gives the other well-known identity 



8. Some arithmetic properties of the sequences a„ and 6„ 

Given a sequence of integers {xn\ it is often interesting to examine its arithmetic 
properties. For instance, given a prime p, this is measured by the p-adic valuation 
I'pixn), defined as the largest power of p that divides a;„. Examples of this process 
appear in [2] for the Stirling numbers and in [1] [22] for a sequence of coefficients 
arising from a definite integral. 

The statements described below give information about Upifln)- These results 
will be presented in a future publication. M. Lasalle [TH] established the next 
theorem by showing that A„ and C„ have the same parity. The fact that the 
Catalan numbers are odd if and only if n = 2'' — 1 for some r > 2 provides the 
proof. This result appears in [T31 [IB] ■ 

Theorem 8.1. The integer a„ is odd if and only if rt = 2(2™ — 1). 

The previous statement may be expressed in terms of the sequence of binary 
digits of n. 

Experimental Fact 8.2. Let B{n) be the binary digits of n and denote x a 
sequence of a arbitrary length consisting of the repetitions of the symbol x. The 
following statements hold (experimentally) 

1) V2(an) = if and only if B{n) = {1, 0}. 

2) V2(an) = 1 if and only if B{n) = {1} or {1,0}. 

3) V2{an) = 2 if and only li B{n) = {1,0,1,0}. 

The experimental findings for the prime p — Z are described next. 
Experimental Fact 8.3. Suppose n is not of the form 3™ — 1. Then 

(8.1) Vi{a'in-2) = I^3(a3ri-l) = t^3(a3ri)- 

Define Wj — ~ 1. Suppose n lies in the interval Wj + \<n< Wj+i — 1. Then 

(8.2) 1^3(0371+2) = j - 1^3 (?^ + !)■ 

If n — Wj, then 1^3 (a3„) — 0. 

Now assume that n = 3™ — 1. Then 

(8.3) l^3(a3n) = i^3(a3n-l) - 1 = J^3(a3n-2) - 1 = m. 



(7.39) 




(7.40) 
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Experimental Fact 8.4. The last observation deals with the sequence {a„(/i)}. 
Consider it now as defined by the recurrence ()7.10p . The initial condition ai(/i) — 2, 
motivated by the origin of the sequence, in general does not provide integer entries. 
For example, if /i = 2, the sequence is 

r 2 8 7 88 1594 14481 
\ ' 3' 9' 3' Y' ~2T~' 

and for /i = 3 

r 1 1 39 263 3091 

\ ' 2' 2' 40' ' ~J ■ 
Observe that the denominators of the sequence for = 2 are always powers of 3, 
but for = 3 the arithmetic nature of the denominators is harder to predict. On 
the other hand if in the case fj, = 3 the initial condition is replaced by ai(3) = 4, 
then the resulting sequence has denominators that are powers of 5. This motivates 
the next definition. 

Definition 8.5. Let Xn be a sequence of rational numbers and p be a prime. The 
sequence is called p-integral if the denominator of a;„ is a power of p. 

Therefore if ai(3) = 4, then the sequence a„(3) is 5-integral. The same phenom- 
ena appeas for other values of ^, the data is summarized in the next table. 





2 


3 


4 


5 


6 


7 


8 




2 


4 


10 


12 


84 


264 


990 


P 


3 


5 


7 


7 


11 


11 


13 



Note 8.6. The sequence {2, 4, 10, 12, 84, 264, 990} does not appear in Sloane's 
sequences list OEIS. 

This suggests the next conjecture. 

Conjecture 8.7. Let /i G N. Then there exists an initial condition ai(/i) and a 
prime p such that the sequence anifJ-) is p-integral. 

Some elementary arithmetical properties of a„ are discussed next. A classical 
result of E. Lucas states that a prime p divides the binomial coefficient (^) if and 
only if at least one of the base p digits of b is greater than the corresponding digit 
of a. 

Proposition 8.8. Assume n is odd. Then a„ is even. 
Proof. Let n — 2m + 1. The recurrence (|2.10p gives 

2m 

2(2m + l)a2m+i = ^ 

For k in the range 1 < fc < m, one of the indices k or 2m + 1 — fc is odd. The 
induction argument shows that for each such k, either Uk or a2m+i-k is an even 
integer. This completes the argument. □ 

Lemma 8.9. Assume n = 2'" — 1. Then ia„ is an odd integer. 



2m + 1\ /2to+ 1 
fc-1 )\ k+l 
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Proof. Proposition 18.81 shows that ia„ is an integer. The relation (jl.8p may be 
written as 

(8.4) (-!)"-«„ ^ 2 + + ' 



This impHes 



(8-5) = ^E(-irG!i)(- + ! 

Observe that if j is odd, then aj is even and ("^x) ^'^ even. Therefore the 
corresponding term in the sum is divisible by 4. If j is even, then Lucas's theorem 
shows that 4 divides ("^|) ■ It follows that the right hand side is an even number. 
This implies that ia„ is odd, as claimed. □ 



The next statement, which provides the easier part of Theorem 18.11 describes 
the indices that produce odd values of a-n- 

Theorem 8.10. If n = 2(2™ - 1), then a„ is odd. 

Proof. Isolate the term j = n/2 in the identity (|8.4p to produce 

|(-ira„ + 21(2"'-l) = (2".-2j( 2" j ^""^ 

- i !:(-)<, !,)(":; 

Lemma 18.91 shows that ^a„/2 is odd and the binomial coefficients on the first term 
of the right-hand side are also odd by Lucas' theorem. Each term of the sum is 
even because aj is even if j is odd and for j even (j"^) even. Therefore the entire 
right-hand side is even which forces a„ to be odd. □ 

The final result discussed here deals with the parity of the sequence &„. The 
main tool is the recurrence 

(8.6) ^"-Ef":')fr;)^^^»- 



k=l 



with 6i = 1. Observe that the binomial coefficients appearing in this recurrence 
are related to the Narayana numbers iV(n, k) (|1.2p by 

(8^7) ("-')(»:;).(„-i)iv(„-i,.-i,. 

Arithmetic properties of the Narayana numbers have been discussed by M. Bona 
and B. Sagan [5]. It is established that if n = 2™ — 1 then N{n,k) is odd for 
< A: < n — 1; while if n = 2™ then N{n, k) is even for 1 < fc < n — 2. 

The next theorem is the analog of M. Lasalle's result for the sequence &„. 

Theorem 8.11. The coefficient 6„ is an odd integer if and only if n = 2™, for 
some m > 0. 



A SEQUENCE RELATED TO NARAYANA POLYNOMIALS 



25 



Proof. The first few terms bi — 1, b2 — I, = A support the base case of an 
inductive proof. 

If n is odd, then 

n-l 

(8.8) bn^{n-l)Y,N{n-l,k-l)bkbn-k 

fc=i 

shows that &„ is even. 

Consider now the case n = 2™. Then Lucas' theorem shows that /r^)(\j]^^) is 
odd for all k. The inductive step states that bk is even if fc 2''. In the case fc = 2'', 
then bn-k is odd if and only if fc = 2™"^, in which case all the terms in (18.81) are 
even with the single expection (j™-!^) (2^-7^1)^2'"-! • This shows that 6„ is odd. 
Finally, if n = 2j is even with j 7^ 2^, then 

^^(^:0(r.>3-£("r)(:::)-- 

Now simply observe that j 7^ 2*", therefore bj is even by induction. It follows that 
bn itself is even. 

This completes the proof. □ 

9. One final question 

Sequences of combinatorial origin often turn out to be unimodal or logconcave. 
Recall that a sequence {xj : 1 < j < n} is called unimodal if there is an index 
such that xi < 2:2 < ■ • • < and Xm,+i > Xm,+2 ^ ■ ■ • ^ The sequence 
is called logconcave if Xn+iXn^i > x'^. An elementary argument shows that a 
logconcave sequence is always unimodal. The reader will find in [H IHl [Zl HI ISl UHl 
[531 US a variety of examples of these type of sequences. 

Conjecture 9.1. The sequences {a„} and {6„} are logconcave. 
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